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Abstract —The merit of evolutionary algorithms (EA) to solve 
convex optimization problems is widely acknowledged. In this 
paper, a genetic algorithm (GA) optimization based waveform 
design framework is used to improve the features of radar 
pulses relying on the orthogonal frequency division multiplexing 
(OFDM) structure. Our optimization techniques focus on finding 
optimal phase code sequences for the OFDM signal. Several 
optimality criteria are used since we consider two different radar 
processing solutions which call either for single or multiple- 
objective optimizations. When minimization of the so-called 
peak-to-mean envelope power ratio (PMEPR) single-objective is 
tackled, we compare our findings with existing methods and 
emphasize on the merit of our approach. In the scope of the 
two-objective optimization, we first address PMEPR and peak- 
to-sidelobe level ratio (PSLR) and show that our approach 
based on the non-dominated sorting genetic algorithm-II (NSGA- 
II) provides design solutions with noticeable improvements as 
opposed to random sets of phase codes. We then look at another 
case of interest where the objective functions are two measures 
of the sidelobe level, namely PSLR and the integrated-sidelobe 
level ratio (ISLR) and propose to modify the NSGA-II to include 
a constrain on the PMEPR instead. In the last part, we illustrate 
via a case study how our encoding solution makes it possible to 
minimize the single objective PMEPR while enabling a target 
detection enhancement strategy, when the SNR metric would be 
chosen for the detection framework. 

Index Terms —OFDM Radar, Genetic Algorithm, NSGA-II, 
PSLR, ISLR, PMEPR 


1. Introduction 

With the advent of powerful digital hardware, software 
defined radio and radar have become an active area of research 
and development [?], [?]. This in turn has given rise to 
many new research directions in the radar community, which 
were previously not comprehensible. One such direction is 
the recently investigated OFDM radar [?] which uses OFDM 
waveforms instead of the classic linear frequency modulated 
(LFM) waveforms. 

In communications, OFDM is a special form of multicarrier 
modulation (MCM), where a single data stream is transmitted 
over a number of lower rate subcarriers. Alphabets such as 
binary phase shift keying (BPSK), quadrature phase shift key¬ 
ing (QPSK), etc., are commonly used to code the information. 
The resulting complex symbols, also called phase codes (- 
l,-i-l in BPSK) are modulating the subcarriers [?]. When the 
signal gets to the receiver, a demodulation stage retrieves the 
transmitted phase codes and eventually the binary message. 
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In radar, the priority is to detect the presence of targets and 
possibly estimate some of their features through the following 
measurable quantities: range, Doppler, azimuth and elevation. 
It is thus of utmost interest for the radar designer to understand 
the possibilities offered by the OFDM structure and optimize 
it to suit its needs. 

It can be quickly verified that non-coded OFDM pulses (all 
phase codes are 1 or at least equal) will not be suitable to 
radar systems that operate with the conventional matched filter 
processing, since they will give rise to high sidelobes. Another 
drawback of non-coded OFDM pulses lies in the behaviour of 
the time domain signal. Strong variations are detrimental since 
they can cause distortions in power amplifiers. Therefore, the 
OFDM pulse needs to be tailored before it becomes a suitable 
radar waveform. 

In this work, we argue that some simple and easy-to- 
implement memetic computing techniques can be used to 
design OFDM radar pulses. Note that the investigation of 
new evolutionary techniques is out of the scope. We look 
at the problem of reducing the peak-to-mean envelope power 
ratio (PMEPR) solely and jointly with the problem of reducing 
the autocorrelation sidelobes. These aspects respond to two 
processing mechanisms specific to OFDM radar that were 
presented in [?]. One of these mechanisms, so called frequency 
domain processing, mostly applies to tracking scenarios. When 
dealing with the sidelobes’ reduction (in regard of the other 
mechanism), because our approach is restricted to the range 
dimension it remains applicable to tracking scenarios where 
the Doppler effect can be accounted for. In this context, one 
would not consider the autocorrelation function but rather, the 
correlation of the pulse with the same pulse including Doppler. 
Although a number of techniques have addressed these objec¬ 
tives separately and jointly [?], [?] they also give rise to several 
limitations. For example, the Newman phasing technique [?] 
gives a PMEPR below 2 for the single OFDM symbol case, 
and this for most numbers of subcarriers N [?, chap. 11], 
but it degrades when the subcarriers are not contiguous as 
we shall see. Another drawback is the impracticality to reuse 
existing communication schemes that, for instance would rely 
on predefined alphabets such as BPSK, QPSK, etc.. The third 
limitation regards the lack of control over the spectrum. This 
is detrimental in the context of matched illumination, which 
is a strong incentive to use OFDM radar [?]. Hence, this 
paper aims to feature emerging evolutionary algorithms (EA) 
in response to the above challenges and includes multiple- 
objective optimization (MOO) considerations, which have not 
been researched before. 

When addressing the single-objective optimization (e.g. 
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PMEPR), we implement the genetic algorithm (GA) in its 
binary [?] or continuous form [?]. The binary form is well 
suited when working with discrete sets of phase codes (i.e. 
constellations). In the case of multiple-objective optimizations, 
we use the well known non-dominated sorting genetic algo¬ 
rithm II (NSGA-II)[?]. It has proven to work much faster than 
the earlier version [?], NSGA, while providing diversity in the 
solutions. Other radar related studies have implemented the 
same algorithm [?], [?], [?]. 

In the context of OFDM signalling for radar problems, 
our contribution lies in the use of genetic algorithm based 
methods, where the parameters are the real arguments of the 
complex phase codes of the OFDM signal, to control the 
OFDM pulse in terms of its time domain variations and its 
autocorrelation sidelobes. To the authors’ best knowledge, 
these considerations are novel and contribute to the field of 
OFDM radar whose state-of-the art as reviewed in [?], [?] is 
still applicable. To compare with more recent works, it must be 
noted that in [?], Sen et al. address a detection problem in the 
condition of multipath, in which MOO-GA based techniques 
are employed to optimally design the spectral parameters of 
the OFDM waveform for the next coherent processing interval. 
There, the objective functions are driven by optimal detection 
considerations and therefore are distinct from those considered 
in this paper. 

The rest of the paper is organized as follows. In Section [n| 
we aim at dimensioning the OFDM pulse that we will then 
optimize, in particular, we assess realistic values for the num¬ 
ber of subcarriers. For this, we briefly review the processing 
solutions proposed in OFDM radar and connect the OFDM 
pulse parameters to external parameters, after formulating 
constraints driven by the processing and the scenario. In 
Section |I^ we review the objective functions at stake in the 
specific case of discrete signals, recall their relevance with 
respect to the processing technique selected and discuss the 
variables involved in the optimization. In Section |IV| we 
present our optimization techniques, starting with a review of 
state-of-the-art methods followed by a discussion on the incen¬ 
tives to employ genetic algorithm based procedures that builds 
upon the limitation of the existing techniques. We then review 
swiftly our encoding strategy. In Section |Vj we present our 
numerical results for the various cases of single objective and 
multiple objectives and establish the merit of the techniques 
by comparing with prior studies, when possible. In Section [Vl| 
we address the specific case of PMEPR optimization as part 
of a matched illumination procedure. We demonstrate that, 
on the one hand, the target detection enhancement problem 
and, on the other hand, the PMEPR optimization problem, 
can be solved if we reserve the real weights to address the 
first problem and use the complex phase codes for the other. 
We give our conclusions in Section |VII| 

II. Dimensioning the OFDM pulse 

In this section, we present the successive steps that form our 
design strategy for fixing some of the OFDM pulse parameters. 
We show how they are inferred from, on the one hand, the type 
of processing and on the other hand, the scenario. This step 


is particularly relevant to reveal some specificities of pulsed 
OFDM radar, as opposed to continuous OFDM radar. The 
results from this section are used later on to generate pulses 
with appropriate dimensions. 

A. Processing related constraints 

In the scope of a pulsed OFDM radar waveform, we 
proposed in [?] two processing alternatives. The first alter¬ 
native is based on the combination of matched filtering and 
Doppler processing while the second alternative transforms 
the received signal in the frequency domain in the same way 
OFDM communication systems operate. After a demodulation 
stage that suppresses the phase codes, two orthogonal discrete 
Fourier transforms (DFT) are applied to form a range Doppler 
image. The key characteristics of both processing are recalled 
in Table |T| In the rest of the paper we refer to the former as 
the time domain processing while we define the latter as our 
frequency domain processing. 

TABLE I 

Characteristics of our OFDM processing alternatives 



Time domain processing 

Frequency domain 
processing 

Pros 

Immune to intersymbol 

Range and Doppler 

interference 

Doppler sidelobes are 
phase codes independent 

sidelobes are phase 
codes independent 

Cons 

Range sidelobes are phase 

Subject to intersym¬ 


codes dependent 

bol interference 


The frequency domain processing is subject to inter-symbol 
interference (ISI). This issue has been of utmost interest in the 
early years of OFDM signalling for communication, to cope 
with the multipath effect. To that end the concept of cyclic 
prefix has been introduced [?]. In pulsed OFDM radar, cyclic 
prefix is no longer a relevant feature, as it is in its continuous 
counterpart [?], hence in this discussion we propose to mitigate 
the risk of ISI by matching the size of the range cell according 
to the target extent. This supposes prior information over the 
target dimension, thus we recommend using this processing for 
tracking. Since the time domain processing does not come up 
with a severe design constraint we choose to use the constraint 
resulting from the frequency domain processing as the main 
guideline to dimension the OFDM pulse. 

1) Sampling frequency: Assuming that we work with the 
complex baseband OFDM signal in the receiver, with band¬ 
width B, the Nyquist theorem states that the sampling fre¬ 
quency can be taken as low as /^ = 5 and the time cell size 
is thus inversely proportional to the bandwidth, tg = l/B. The 
size of the range cell is then given by c/2B, where c is the 
speed of light. 

2 ) Bandwidth: To satisfy the above constraint on the size of 
the range cell we shall adjust the bandwidth to comply with: 
c/25 > A5t where A5t is the target range extent. Practically 
we can add a margin to account for the target radial velocity 
and the uncertainty on the target extent and position. Figure 
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illustrates this constraint, which is itself inferred by the loss of 
orthogonality of the subcarriers caused by ISI [?]. Not only the 
received echoes shall remain in the same time cell as a result 
of the first pulse but also throughout the coherent processing 
interval. In the end the bandwidth can be obtained from: 


B = 


■ a) 


( 1 ) 


Fa\ ourable case 
__ C/2B .— 


2(Ai?, 

where a is the margin in range. 

Unfa\ ourable case 
_. C/2B . _ 

Echo scatterer #1 

Echo scatterer #2 

Echo scatterer #3 


Fig. 1. Constraint applicable to the size of the range cell in order to mitigate 
the risk of ISI in the case of the frequency domain processing. Echoes received 
from the different scatterers that compose the target shall fall within the same 
range cell. 




to the bandwidth and the maximum number of subcarriers for 
both cases. 

TABLE II 

SCENARII CHARACTERISTICS EOR THE WAVEEORM DESIGN 



Case 1 

Case 2 


(walker) 

(truck) 

Range extent (m) 

2 

10 

Margin (m) 

1 

5 

Bandwidth (MHz) 

50 

10 

Maximum number of subcarriers 

500 

100 


III. Optimizing the pulse for radar 

Having dimensioned our OFDM pulses we are now ready to 
concentrate our analysis on the optimization of those OFDM 
pulses for radar. Firstly, we need to define our objective 
functions and secondly, we need to identify the OFDM pa¬ 
rameters or variables upon which we can run our optimization 
procedure. 


Therefore, if there is no need for high range resolution 
(assuming that the target is known) we suggest to dimension 
the bandwidth based on these ISI considerations. 


B. Scenarios related constraints 

Other parameters need to be fixed, viz., the pulse length and 
the number of subcarriers. 

1) Pulse length: To dimension the pulse length, we relate 
it to the eclipsed zone, the window that lies between the radar 
and the minimum detection range like it is connnonly done for 
pulsed radar. If, say, we expect targets from = 1.5 km, 
an upper bound for the pulse length is found to be [?], tp = 
2-Rmin/— 10 

2) Number of subcarriers: The orthogonality property is 
another example of the unique OFDM structure. It states that 
the symbol or bit duration 4 is inversely proportional to the 
subcarrier spacing A/, R = 1/Af. In the extreme case where 
the pulse is composed of one symbol, the maximum number 
of subcarriers A^ax used in the pulse can be derived from: 

2Bi?niin 


Aimax - 


( 2 ) 


A. Objective functions 

An OFDM symbol is built as a sum of weighted complex 
sinusoids, where every sinusoid has a given starting phase. Be¬ 
cause the multiple frequencies (or subcarriers) are transmitted 
simultaneously the full bandwidth is covered at every instant 
in the OFDM configuration [?]. Fig. [^depicts an OFDM pulse 
and compares it with the typical LFM pulse. 




where we used the previous result for the pulse length upper 
bound. For the same pulse bandwidth and pulse duration a 
smaller number of subcarriers can be used if we construct the 
pulse from several symbols. For example, we can decide to 
use 250 subcarriers with 4 symbols in the pulse to maintain 
the same duration. The subcarrier spacing is then increased 
from 100 kHz to 400 kHz. Information related to the target 
maximum speed shall also be considered to ensure negligible 
distortion effects in either processing [?], [?]. 

To put the above analysis into perspective, we suggest two 
objects of different size, on the one hand a walker and on the 
other hand, a truck. Table |n| summarizes the values attributed 


(b) OFDM pulse 


Fig. 2. Real part of a) a LFM pulse versus b) an OFDM pulse with random 
starting phases. Both pulses have been normalized to have unit energy. 


When the OFDM pulse is composed of several symbols it 
can be expressed as: 

N / K \ 

Wn i'^an,krk{t) \ exp{j2TTnAft). (3) 


n=l 


Kk=l 


A is the total number of subcarriers in the pulse, denotes 
the weight applied on subcarrier n and an,k is the phase 
code attributed to subcarrier n in symbol /c. A is the total 
number of symbols in the pulse. The function rk{t) refers to 
the rectangular window for every symbol: 


rk{t) = 


{i 


kR <t < {k ^ 1)4 

elsewhere 
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Like in [?], the normalization factor A is given by: 

1 


y^K-l y^N-1 2 I 

Z^/c=0 Z^n=0 ^nl 


^n,k 


\Ht 


( 4 ) 


When using OFDM for radar two important aspects must 
be considered. 

1) Peak-to-mean envelope power ratio: This aspect is spe¬ 
cific to OFDM signals. It characterizes the variations in time 
of the signal envelope as observed in Fig. |2(b) Because they 
are generally used in saturation to maximize the transmitted 
power, devices like power amplifiers can distort the signal 
and be detrimental to the system performance. While the 
peak-to-average power ratio (PAPR), defined as the ratio of 
the peak power to the average power of the real-valued 
multicarrier signal, is usually of interest, we consider the peak- 
to-mean envelope power ratio (PMEPR) because we work with 
complex baseband signals. By definition, the PAPR will result 
in smaller values than the PMEPR [?] whose expression is 
given by [?]: 


max X n 


PMEPR = 




( 5 ) 


2) Sidelobe level: Typically, with the time domain pro¬ 
cessing where matched filtering is applied in range, it is 
important to keep the sidelobes of the compressed pulse as 
low as possible. Hence, our next objective functions regard the 
behaviour of the autocorrelation function, whose expression 
R{t) is given by: 


R{t) 


f 


— r)dt. 


and in its discrete form, R[r 


R[m] = 


NK-l 

E 

p=0 


x[p]x*[p — m] 


( 6 ) 


( 7 ) 


where m takes integer values between —NK +1 and NK — 1. 
X [p] represent the discrete values of the OFDM pulse taken at 
the discrete instants pR/N, where p takes integer values from 
0 to NK — 1. If the pulse is composed of only one symbol, 
then p takes values from 0 to A" — 1, just like n the subcarrier 
index. In the end, x[p] = x{ptb/N). In Eq. we assume that 
x[p] = 0 for all forbidden values of p, that is p < 0 and 
p > NK — 1. 

To cope with practical applications we commonly distin¬ 
guish two objective functions. The first function is the PSLR. 
It returns the ratio between the highest sidelobe and the peak. 


max I i?[m] | 

PSLR= 


|i?[0]| 


( 8 ) 


The second function is the integrated-sidelobe level ratio 
(ISLR). It returns the ratio between the cumulation of the 
sidelobes and the peak. 

EI-RNI 


ISLR = 


|i?[0]| 


-,m 7^ 0 


( 9 ) 


The weight attributed to both figure-of-merits depends on 


the application as well as the environment. For example, if 
the radar operates in the presence of distributed clutter, it 
will be important to work with low ISLR in order to keep 
the weak targets visible. In that case, high ISLR can be 
interpreted as an increase of the noise fioor. Conversely if 
the application requires detection of targets in the presence of 
strong discrete clutter, the PSLR is more critical and must 
be kept low to prevent from deceptively considering one 
sidelobe as another small target. In Section [V| we present our 
results using the conventions PSLR^^ = 20\og{PSLR) and 
ISLR^b = 20log{ISLR). Remark that the above objective 
functions are for the zero Doppler case. In fact, would the 
design concern a tracking scenario, the Doppler effect could 
be accounted for in the model of the received signal such that 
the objective function would no longer be inferred from the 
autocorrelation function but simply the correlation of the pulse 
with the pulse including Doppler. In the case where there is no 
a priori information concerning Doppler, the optimization shall 
be done over the range Doppler plane that includes all possible 
Doppler values. In [?] we showed that a 1 dB compression loss 
on the main peak would arises when jNf /K, hence if the 
subcarrier spacing is large enough (few MHz) the OFDM pulse 
is Doppler tolerant. However, this does not conclude on the 
sidelobe level over the range Doppler domain. Considerations 
regarding the optimization of the range Doppler sidelobes are 
outside the scope of this paper. 

3) Need for oversampling: In Eq. 0Eq.@ and Eq. we 
have assumed to work at the critical sampling rate fs = B 
such that the sampling period is tg = h/N, as a result of 
the relationships that govern the OEDM structure. In Eig. 
we emphasize on the need for oversampling. Because of the 
quick temporal variations of the OEDM signal, the PMEPR of 
the critically sampled signal will differ from the PMEPR of 
the continuous signal. An important question is how large the 
oversampling factor should be in order for the approximation 
to be accurate. In [?], the difference between the continuous 
time and discrete time PAPR is evaluated and the conclusion 
is that an oversampling factor of 4 is accurate for an OEDM 
signal with BPSK coding. In this work, we choose to oversam¬ 
ple the OEDM signal with a factor of 20 before evaluating our 
objective functions. In Eig. |3(b) we show how both PSLR and 
ISLR calculations exclude the values around the main peak. 
The total peak extent is equal to twice the Rayleigh resolution, 
which in time is 2/B [?]. 

B. Single or multiple objectives 

In light of the above, we see that in the process of designing 
an OEDM pulse, extra care shall be employed to maintain tol¬ 
erable PMEPR as well as reasonable sidelobe levels. However, 
if we rely on our frequency domain processing, we have re¬ 
called in Table 1^ that both the range and Doppler sidelobes [?] 
were phase code independent. In this case, the pulse design can 
focus on the PMEPR only. In a nutshell, with pulsed OEDM 
radar, we suggest two distinct strategies for the design of the 
pulses, viz., a single objective PMEPR optimization for the 
frequency domain processing and a multi-objective PMEPR 
and autocorrelation sidelobes’ level optimization for the time 
domain processing. 
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(a) Time domain OFDM pulse (b) Autocorrelation output 

Fig. 3. a) Amplitude versus time of an OFDM pulse and b) its autocorrelation 
function versus range when N=3 and K=3 in the case of critical sampling 
fs=B and oversampling fs=20B. The values obtained for the PMEPR and 
the PSLR are given. 

C Variables for the optimization 

The parameters that remain after dimensioning the OFDM 
pulse are, as seen in Eq[^ the weights Wn and the phase 
codes an,/c. Together, they form a vector of complex numbers 
u;nexp(j0n,/c), where an,/c=exp(j0n,/c) and fn.k are real 
numbers between 0 and 2tt. As mentioned earlier, the weights 
can also be regarded as spectral coefficients and as a result of 
the scattering mechanism observed in radar, their value may 
infiuence the detection process. Hence, it may be favourable to 
disregard them as variables for the optimization of the above 
objective functions. 

IV. Optimization techniques 

In this section, we first revisit state-of-the-art techniques 
that addressed the objective functions presented earlier. Then, 
we discuss their limitations and motivate the use of genetic 
algorithm based methods. 

A. State-of-the-art methods 

The problems of minimizing the PMEPR and the PSLR 
have been extensively researched since the emergence of 
multicarrier radar signals. Two strategies have stood out [?, 
chap. 11]. The first approach is to consider identical sequences 
(IS), such that all subcarriers are assigned the same phase code. 
Optimizing the PMEPR of the pulse results in optimizing the 
PMEPR of a single symbol. In this context, Newman [?], 
Schroeder [?] and Narahashi [?] have suggested different 
phasing methods, where the weights Wn are complex values 
with quadratic dependence on n. PMEPR as low as 2 (3 dB) 
can be obtained for any relevant value of N (up to 65,000). 
Sidelobe levels of the autocorrelation function (ACF) can 
then be adjusted by introducing a linear carrier phase term 
that does not impinge the PMEPR but can change the ACF 
sidelobe pattern [?]. The second design strategy is based 
on modulating all N subcarriers with consecutive ordered 
cyclic shifts (COCS) of an ideal chirplike sequence (CLS) 
of length AT, assuming that K symbols compose the pulse. 
For example, an OFDM pulse based on COCS of P4 codes 
can give PMEPR below 2 and PSLR below -15 dB for a large 
range of N (between 0 and 70). Besides, note that a handful of 
methods aiming to reduce the PMEPR of OEDM symbols for 
communication have been developed [?], however they refer 


to the transmission of information and hence fall outside the 
scope of this discussion. 

B. Motivation for genetic algorithm optimization methods 

Although the aforementioned techniques produce excellent 
results in terms of PMEPR alone or PMEPR and PSLR 
combined, they suffer from several limitations that may be 
detrimental in prospective situations of interest. Eirstly, those 
techniques need contiguous subcarriers: from 1 to N. Hence, 
what if we want to use a pulse with a sparse spectrum 
(as a result of jamming, interference or simply frequency 
allocation) [?]? The objective functions will be degraded un¬ 
doubtedly. Secondly, they imply that existing communication 
schemes (such as BPSK, QPSK etc. constellations) cannot 
be used, which may prevent from building upon existing 
communication infrastructure to perform radar. Thirdly, the 
weight vector is fixed in the optimization mechanism of those 
objective functions. This prevents from shaping the spectral 
components of the signal to improve the detectability of an 
object whose radar cross section (RCS), we know, resonates 
variably at different frequencies [?, chap. 7], [?]. In this 
context, GA based methods offer favourable optimization 
strategies for they are able to address these aspects and are 
easy to implement. Other incentives include the possibility 
to work with modified objective functions, add new objective 
functions such as ISLR, propose a multitude of ’good’ solu¬ 
tions and provide Pareto fronts of optimal solutions in a multi¬ 
objective sense. Until now, just a few authors have considered 
the use of GA and MOO techniques for radar, however an 
increasing number of papers have recently established the 
potential of these optimization methods, in the context of 
improved detection under multipath conditions for OEDM 
radar [?], range ambiguity suppression in OEDM SAR [?] and 
OEDM-STAP detection [?]. Other optimization techniques that 
are receiving increasing attention within the broad engineering 
field are convex methods [?]. Although we do not claim 
that our problem is not convex or cannot be made convex, 
considerations regarding the use of convex optimization is 
beyond the scope of this paper. 

C. Problem encoding 

In our approach, the variables consist of the NK arguments 
(real numbers) Note that the OEDM symbols result from 
applying an inverse discrete Eourier transform (IDET) on each 
of the N long vectors [wn G^'p{jfn,k)]n=i---N [?]• The OEDM 
pulse is then obtained by concatenating the symbols if AT > 1. 
The variables can either be encoded in the form of binary 
strings [?] or alternatively real values [?], [?]. In the case 
where we would employ a predefined constellation, it makes 
sense to use the first option since ’by nature’, each phase 
code translates into a short string of few bits (e.g. two bits 
for QPSK, 3 bits for 8-PSK etc.). In the continuous case, 
the resolution depends on the system resolution, which in our 
case, allows us to represent one real number with 18 bits. The 
resolution in angle is then AO = 27r/2^^ 0.024 mrad. Note 

that we choose 18 to maximize the size of our domain of 
solutions and hence increase the chance to approach the best 
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solution at the cost of increasing the search space. In practice, 
one may want to restrict it to a smaller number for the sake 
of reducing the computation complexity. 

V. Numerical results 

In this section we present our results for the various config¬ 
urations, single and multiple-objective optimizations. We start 
with the case of the single objective PMEPR. Our intention 
is to show that the genetic algorithm is an elegant method to 
find optimal sets of phase codes, while challenging important 
limitations of state-of-the-art techniques. We then present our 
results for the case of multiple objectives and develop our 
discussion to introduce a PMEPR constrained multi-objective 
optimization where the original NSGA-II is slightly modified. 

A. Single objective: PMEPR 

In our simulations, we implemented a simple genetic algo¬ 
rithm (SGA) [?]. We considered a population of 12 chromo¬ 
somes [?] and included mutation (Hipping one bit of the binary 
string that forms the chromosome, at random) at the rate of one 
mutation for 5 ’offspring’ chromosomes, selected at random 
every two generations. Despite the simplicity of our genetic 
algorithm, we are able to retrieve phase code sequences with 
optimal PMEPR properties. 
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Eollowing, the preliminary analysis from Section N=500 
can be considered as an upper bound for pulsed OEDM radar. 
Hence, most likely, smaller number of subcarriers may be 
used, inferring computation times less than 500s. Whether the 
search can be implemented on the fiy or not will depend on 
the infrastructure available in terms of the computation power. 
Otherwise, it may be done off line. 

Pigs. 5(a)|5(c) illustrate the complex phase codes corre¬ 
sponding to one chromosome of the final population for 
both cases. It is clear that a pseudo-continuum of values are 
available for the first case, while with QPSK, the arguments 
(j)n,k are limited to 4 values: 0, 7r/2, tt and 37r/2. 
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Fig. 4. Convergence of the SGA while encoding one phase with either a) 
18 or b) 2 genes, for N=100 and N=500. A minimum PMEPR of 2.9 is 
achieved after 400 generations in average for both cases, when N=100. In 
fact, we identified a PMEPR as low as 2.6 and 2.7 in our populations after 
400 generations, respectively for the M-PSK (M=2^®) and QPSK cases. 

We have run a Monte-Carlo simulation with 20 runs and for 
each run at each generation, we found the minimum PMEPR 


in the population. The values reported in Pig. |4(a)| result from 
averaging out these minima over the multiple runs. This means 
that, for the case where A'=100, after 400 generations in 
average one chromosome of the final population (one set of 
phase codes) generates an OEDM symbol whose PMEPR is 


2.9. While Pig. 4(a) shows the c onver gence of our GA for 
the case of M-PSK (M=2^^) Pig 4(b) refers to the case of 
QPSK modulation. With Pig |4(b)[ we demonstrate that it is 
possible to control the PMEPR of an OEDM symbol that 
would rely on a predefined alphabet. In terms of the execution 
time, we found that in average, it took 115s to loop over 
400 generations for the first case, when it took in average 
104s for the second case, with A^=100. When A^=500, this 
time increases to 566s and 515s respectively. We are running 
our simulation on the university of Cape Town (UCT) server. 


Fig. 5. Illustrating the complex phase codes and sparse spectrum in the 
case of, a) and b) M-PSK (M=2^®) and sparsity coefficient 50% and in the 
case of, c) and d) QPSK and sparsity coefficient 70%. We consider N=100 
subcarriers. 

The second limitation that we challenge with our GA based 
method relates to the case of sparse spectra. To that end, we 
repeated the above simulations for sparsely populated spectra, 
with sparsity coefficients of 70% (out of 100 subcarriers, 70 
subcarriers are ON, 30 are OPE) and 50%. At each run of 
the Monte-Carlo simulation, a sparse spectrum is selected at 
random. In practice, we populate the band at random but 
maintain the two extreme subcarriers ’ON’ to prevent from 
reducing ’badly’ the bandwidth. Examples of spectra are given 
in Pigs. |5(b)|5(d)| for the two sparsity options. Remark that, 
on the one hand, the weighting is uniform for all subcarriers 
’ON’ and on the other hand, its level is adjusted to provide 
an OEDM symbol of unit energy. 

Por comparison purposes, we assessed the behaviour of 
the Newman phasing technique under the same conditions 
of sparse spectra. We also ran Monte-Carlo simulations and 
collected the mean value for the PMEPR over all runs. Taking 
the mean value makes sense since the decision as to which 
subcarriers shall be OPE may not be in the designer hands. 
We considered 1000 runs. Our results, given in linear unit, are 


presented in Table III for both cases A^=100 and A^=500. We 
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observe that, under sparsity conditions, they outperform, in 
average, those obtained with Newman’s technique. Note that, 
with A^=500, we achieved better results in the constrained case 
of QPSK as opposed to the ’no constellation’ case, after 400 
generations. One cause is the size of the search space, much 
bigger in the unconstrained case. Whether we can improve 
these results further would require optimizing our genetic 
algorithms for the case of large N. However, this falls outside 
the scope of this work. 

TABLE III 

Simulation results for the single objective PMEPR (in linear 
unit) with different levels of sparsity. Our GA based results 

ARE COMPARED WITH TWO KNOWN RESULTS, VIZ., THE TRIVIAL 
NON-CODED CASE AND NEWMAN PHASE CODES’ CASE 


We also show that our GA based technique can even outper¬ 
form Newman’s phasing method in the case of fully populated 
spectrum, with smaller number of subcarriers, viz., ^"=10 in 
this example. Fig. |6(a)| presents our PMEPR optimized OFDM 
symbol and compares it with the Newman based symbol. 
Fig. |6(b)| shows another PMEPR optimized OFDM symbol 
in the case of 80% sparsity. The difference between both GA 
based and Newman based PMEPR is about 1 dB. 
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that the PMEPR optimization does not impinge the matched 
illumination and reciprocally if they are addressed in the 
correct order. 

B. Multiple objectives 

One benefit to use the NSGA-II in multi-objective optimiza¬ 
tion procedures lies in the fact that all objectives are optimized 
concurrently. The Pareto front created at each generation 
converges towards some optimum that provides the designer 
with the best set of solutions in a Pareto sense [?]. In our 
OFDM pulse design problem, needs for a multi-objective 
optimization arise when the time domain processing would 
be selected as explained in Section |III-B| 

1) PMEPR and PSLR: We start by looking at the two- 
objective optimization with PMEPR and PSLR. In Figs. |7(a)[ ^ 
7(b) we show how the NSGA-II can improve both objectives. 
We compare the Pareto front after 10000 generations to a 
random population. In both cases, we have taken a num¬ 
ber of subcarriers smaller than 100 and 500, but the time 
bandwidth product remains the same, ^"7^=100 and NK=500 
respectively. The size of the population is equal to 40. After 
10000 generations our set of optimal solutions is considerably 
improved. The Pareto front is visible and different sets of phase 
codes can be selected depending on whether the emphasis in 
on a low PMEPR or a low PSLR. In light of Eig. |7(a)| we see 
that improvements as high as 7 dB and 4 dB can be achieved 
in terms of PSLR and PMEPR for the case of N=25, K=4. 


Likewise, Lig. |7(b)| shows that improvements as high as 9 dB 
and 3 dB can be achieved in terms of PSLR and PMEPR for 
the case of A^=125, K=4. The improvements reported here 
refer to the characteristics of the mean point in the random 
distributions. 


- GA, PMEPR=1.75 

— Newman, PMEPR=1.84 
.No coding, PMEPR=10 




l 0.4 0.6 

Time (//s) 

(a) Full band 



(b) 80% of the band 


Fig. 6. Amplitude of the time domain OFDM symbol when a) all subcarriers 
are used and b) only 80% of the subcarriers are used, N=10. In both cases, our 
GA solution outperforms the Newman strategy. All pulses have unit energy. 


The third limitation raised by state-of-the-art techniques 
regards the lack of control over the spectrum, viz., the weights 
cannot follow any arbitrary pattern. We already observed this 
effect with the degradation of the results following Newman’s 
technique in the case of sparse spectra. Sparsity is, to a certain 
extent, a case of non uniform weight pattern. Although the 
GA based results presented so far have been obtained with 
equal weights for all subcarriers ON, one could easily repeat 
the simulations including non uniform values. We propose to 
address separately this aspect in the context of a matched 
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Fig. 7. Comparison between random phases and phases resulting from the 
NSGA-II optimization in a) and b). Convergence of the NSGA-II optimization 
at constant time bandwidth products in c) and d). 


illumination procedure in Section VI Our aim is to show We then propose to assess the position of the Pareto front on 





















































the PMEPR-PSLR map for different number of subcarriers N 
and compare three cases with equal time bandwidth products. 
We report the results in Figs. |7(c)f |7(d)| For each case, we 
plot the Pareto front after 1000 and 10000 generations. We 
note that designs with small N will provide the best Pareto 
front. In these two-objective simulations, a cycle of 1000 
generations takes about 70s and 370s for the cases, N=25 
K=4 and A^=125, K=4 respectively. 

2) PMEPR constrained multi-objective optimization: One 
may argue that it is more appropriate to treat the PMEPR as a 
constraint rather than an objective. Indeed, if the system pro¬ 
vides a maximum value for the PMEPR, this may leave space 
for other objective functions, like the ISLR. As mentioned 
before, PSLR and ISLR pick up different behaviours of the 
sidelobe level whose impact will be scenario dependent. It is 
thus judicious to generate phase codes that can optimize both 
objectives in the form of Pareto fronts instead of treating each 
objective separately [?]. We report our preliminary scheme to 
deal with this constrained multi-objective problem as follows. 
Our idea has been to modify slightly the original version of the 
NSGA-II and in particular assign very low crowding distances 
to those chromosomes that result in PMEPR exceeding the 
threshold that we chose to apply. The motivation here is to 
reduce the chance for those chromosomes to participate to the 
evolution, such that over time, the population is composed 
of chromosomes with favourable PMEPR values. Although 
quite simple, this approach proved to work on the few cases 
analysed. Hereafter, we present the results of our analysis for 
the configuration A^=100, K=\An order to select a reasonable 
threshold value, we first assessed the PMEPR distribution, 
which we display in Fig. |8(a) We selected the PMEPR 
value preceding the peak, 5 in this case. Our Monte-Carlo 
executed 100 runs of 1000 generations each and we analysed 
the characteristics of all 100 last populations. We considered 


the same parameter for the population size, 40. In Fig. 8(c) 


we assess the PMEPR values for all elements of the last 
populations. We see that a big majority of chromosomes have 
their PMEPR below the threshold. Counting those populations 
that gave rise to 40 chromosomes with a PMEPR below the 
threshold, we get 38 cases out of 100. We can then plot 
the Pareto fronts after 1000 generations for all these 100 
runs, distinguishing those fronts where the constrain is fully 
compliant and the other. We obtain the results as presented 


in Fig. |8(d)| Comparing with the case where the optimization 
is run without constrain on the PMEPR, whose results are 


given in Fig. 8(b) we see that our constrained optimization 
gives Pareto fronts with reduced PSLR values but equally good 
ISLR values. 

For completeness, we repeated our analysis for the case 
A^=500. The threshold for the PMEPR, selected in the same 
manner was chosen to be 6.5. Among the 100 final popu¬ 
lations, only 22 fulfilled the criteria whereby all 40 chro¬ 
mosomes shall induce a PMEPR below 6.5. Comparing the 
results in the unconstrained case in Fig. |9(b)| with those from 
the constrained case as in Fig. 9(a) we see that again, our 
constrained optimization gives Pareto fronts with similarly 
good ISLR values. The performance is not as great for the 
PSLR since the improvement is reduced from 6 dB down to 
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Fig. 8. a) Distribution of the PMEPR for random phase codes, N=100 b) 
100 Pareto fronts obtained after 1000 generations with NSGA-II, no constrain 
on the PMEPR compared with random sets c) distribution of the PMEPR 
for all chromosomes from the last populations after 1000 generations with a 
constrain on the PMEPR d) 100 Pareto fronts obtained after 1000 generations 
with NSGA-II, with a constrain on the PMEPR equal to 5. The Pareto fronts 
where the constrain has been successful for all chromosomes of the population 
are in blue (38 out of 100). 


0.75 dB. 



00 


5 48 
46 
44 


0.75 dB 


dB 

Constrain NOT successful' 

Constrain successful 



-20 -18 -16 -14 -12 

PSLR (dB) 

(b) With PMEPR constrain 


Fig. 9. a) 100 Pareto fronts obtained after 1000 generations with NSGA- 
II, no constrain on the PMEPR compared with random sets, N=500 b) 100 
Pareto fronts obtained after 1000 generations with NSGA-II, with a constrain 
on the PMEPR equal to 6.5. The Pareto fronts where the constrain has been 
successful for all chromosomes of the population are in blue (22 out 100). 


VI. Case study: CA for PMEPR optimization in a 

MATCHED ILLUMINATION PROCEDURE 

Matched illumination is a research direction that motivates 
OFDM radar [?]. In this section, our objective is to show how 
single objective CA based optimizations can be incorporated 
in a matched illumination procedure. We propose a two-step 
approach that addresses both problems of detection enhance¬ 
ment and PMEPR optimization sequentially. As such, this 
problem does not require having a multi-objective approach. 
The first step consists in finding an optimal weight vector 
v^opt to enhance the detection characteristics. In other words. 
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it finds the power spectrum of the radar pulse. This vector is 
then passed on to our single objective GA based optimization 
to find a vector of phase codes aopt (only one symbol is 
considered, hence we drop the index k) such that the resulting 
OFDM symbol has optimal PMEPR properties. Our procedure 


is presented in Fig. 10 The complex vector resulting from the 
element-wise multiplication Wopt • aopt represents the discrete 
spectrum of our optimal OFDM pulse, which we define as 
Xoptif) in Fig. 


10 


X(f) ■ 


enhanced target 
detection 


■ I ■ 


improved PMEPR 
with our GA based 
optimization 




Fig. 10. OFDM pulse optimization in two steps. 


To understand how the weights of the OFDM signal can be 
used to enhance the target detection, we shall mention that in 
radar, scattering can be modelled as a linear filtering process 
where the received echo results from the convolution of the 
transmitted pulse with the target impulse response. Looking 
at it in the Fourier domain, we see that the spectrum of the 
received echo will follow from multiplying the spectrum of the 
transmitted pulse with the target refiectivity spectrum. Hence, 
the shape of the transmitted pulse spectrum, determined by 
the weights will impact the energy content of the received 
pulse. This last inference results from ParsevaTs theorem that 
links the energy of a time domain signal with the energy of 
its Fourier transform in the frequency domain. 

A. SNR as the design metric 

Unlike the traditional approach to consider imaging based 
metrics for target detection and classification. Bell has in¬ 
troduced the signal-to-noise ratio (SNR) metric [?]. In that 
concept, the objective is to maximize the energy or SNR of 
the received echo for the particular target in the scene. The 
PSLR and ISLR are no longer relevant and can be ignored in 
the pulse design. In the following we revisit the formulation 
of this SNR maximization problem. Note that an important 
aspect concerns the definition of the matched filter. 

1) Transmitted signal receiver-filter pair: In [?] the authors 
expressed the signal at the output of the matched filter V{f) as 
the product of the target refiectivity spectrum ^(/), the pulse 
spectrum at baseband X[f) and the filter transfer function 

H{f): 

V{f)=<;if + fc)X{f)H{f). (10) 

In the frequency domain, the pulse goes through the steps as 
indicated in Fig. Both convolution symbols charaterize the 
up and down conversions. In the time domain, we can express 
the radio frequency (RF) analytical signal as 

Vtxit) = x{t)exp{j2TTfct), (11) 

where fc is the carrier frequency and x{t) is the OFDM pulse 
as given in Eq. 

2) Optimum receiver-filter pair: In the case of an extended 
target, unlike the standard matched filtering problem applica¬ 
ble when the target is modelled as a point scatterer, the optimal 


Filter 
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c(/)| \ Sn„(0 

X(f) 5(C ^0 ^ ^ 

t \ 
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Fig. 11. Model of the received signal in the frequency domain. 


receiver filter shall be matched to the waveform scattered by 
the target, not the transmitted target itself. When this is the 
case, the maximum SNR at time to is given by the following 
equation [?]: 

k(/ + /c)^(/)P 


to 


Snnif) 


-df. 


( 12 ) 


While previous works have researched solutions to this SNR 
maximization problem in the time domain, in terms of x{t), 
we propose to run the search in the frequency domain and look 
for |^(/)| instead. This finally amounts to searching for an 
optimal weight vector Wopt- If we assume that the spectrum of 
our OFDM pulse is zero outside the bandwidtlQ and consider 
that the stationary additive Gaussian noise is white, with one¬ 
sided power spectral density Snn{f) = Xq, we can modify 
Eq. p^into: 


to 


TjJ \X{f)c;{f + f,)\^df, (13) 


Assuming further a deterministic expression for the target 
reflectivity spectrum, we can rewrite Eq. f^in discrete terms: 


N-l 


^ to n=0 


(14) 


where X[n]=X{nAf) and <;[n]=g{nAf fc) are the baseband 
OEDM spectrum and the target refiectivity spectrum at RE, 
respectively, expressed in discrete terms. The index n takes 
values from 0 to N — 1. Thus, our optimization problem 
reduces in finding the weights Wn according to: 


N-l 


arg max > u;^k[n]p, s.t. 


n=0 


x'^x= 1 

V\<Wn< Vu 


(15) 


The first constrain tells that the pulse energy remains unity 
throughout the search and the second constrain limits the 
dynamic range of the weight values. Different values for 
the lower and upper bounds will result in different spectrum 
characteristics. Although other methods may exist to solve this 
constrained single objective problem, we suggest applying our 
GA approach, here again. The expression with the summation 
is the new objective function and the variables are no longer 
the phase codes but the weights. 


^The spectrum of the baseband OFDM pulse is considered to spread from 
0 to B rather -BH to BH 
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B. Simulation setup and results 

In this analysis we invoke different dimensions for our 
OFDM pulse as compared to Section [V] We assume to work 
at X-band and consider a 2 GHz bandwidth B centred around 
10 GHz, fc = 9 GHz. Our OFDM pulse is composed of A^=100 
subcarriers and the subcarrier spacing is A/=20 MHz. To 
appreciate the detection enhancement resulting from the op¬ 
timal weight vector, we must normalize the target reflectivity 
spectrum <;[n]. 

1) Target reflectivity spectrum normalization: We normal¬ 
ize the discrete target reflectivity spectrum following the 
strategy presented in [?, chap. 14]. The later states that when 
a flat spectrum OFDM pulse of unit energy interacts with 
the normalized reflectivity spectrum ^norm[^], the frequency- 
domain reflected signal shall have unit average power, viz., 
1/^ • Sn=o^^nknormNP=l. Thc wcights corrcspondii^ to 
our flat spectrum OFDM pulse of unit energy are giveij^by 
Wn = 1/VN. Hence, we And: 

N 

^norniN = . = • g[n], (16) 




2) Complex target model: We simulate a synthetic target 
composed of P=50 point scatterers, each with the same 
reflectivity Q=y/o7=l and located at ranges Ri from the radar. 
The individual point scatterers are assumed to be perfectly 
conducting spheres, large enough to have a reflectivity constant 
within the band of interest. It can be shown [?] that the 
compound target reflectivity spectrum g{f) is equal to: 


Ri. 


?(/) = X 


(17) 


The point scatterers are randomly distributed within a rectan¬ 
gle, 5 meters wide and 10 meters long, whose center is 10 km 
away from the radar along the x axis. 

3) SNR and PMEPR improvements: Following the above 
methodology, we first generate an optimal weight vector Wopt- 
We employ a continuous GA [?] with 20 chromosomes, a 
mutation rate of 0.2 and select arbitrarily i;i=0.01 and ^u=10. 
We also include in the initial population a scaled version of 
the vector <^norm, since we know that this is a good solution. 
We stop the GA after 5000 generations, and consider the final 
weight vector to be Wopt. This takes about 50s. Fig. 12(a) 


table IV The gain is a measure of the difference in average 
power with the case of a flat OFDM spectrum and is equal to 
2.8 dB. This value results from averaging over 10 runs in a 
Monte-Carlo simulation. 

Now that we have fixed |Xopt| the second step consists 
in feeding Wopt into our second GA where the objective 

^In fact, different values are possible depending on how one implements the 
IDFT algorithm, and in particular whether a multiplicative constant is applied 


function is the PMEPR. We consider again 12 chromosomes 
in the population and run it for 600 generations. This takes 
about 170s. A measure of the convergence of our solution is 


presented in Fig. 12(b) We And the best element of the final 
population and define it to be aopt- We observe that the PMEPR 
has been improved by 3 dB. This results is in agreement with 
our results from Section V-A In the end, as shown in Fig 
our pulse Aopt(/) has enhanced detection capabilities for the 
target of interest and optimal PMEPR. 
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(a) Optimal weight vector 
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Fig. 12. a) Optimal weight vector Wopt versus the normalized reflectivity 
spectrum resulting from a complex target made of P= 50 point scatterers and 
b) convergence of the single GA based technique to And an optimal set of 
phase codes aopt- 


TABLE IV 

Result for the detection enhancement. The gain is measured in 

TERMS OF AVERAGE POWER 



Flat spectrum 

GA based Wopt 

Gain (dB) 

0 

2.8 


puts on the same graph the normalized reflectivity spectrum 
for our complex target as well as the optimal weight vector, 
scaled to fall within the same boundaries. We see that most of 
the energy is distributed on those subcarriers that correspond 
to frequencies where the reflectivity spectrum is strong. We 
present the gain resulting from this optimal weight vector in 


VH. Conclusion 

In this paper, we showed that GA based techniques are 
suitable to optimize or improve several features of the OFDM 
pulses, provided reasonable search space, hence, reasonable 
numbers of subcarriers (and symbols). In pulsed OFDM radar 
though, we expect limited numbers of subcarriers as a result of 
the orthogonality principle that links the symbol duration to the 
subcarrier spacing. We inspected the possibility to incorporate 
these optimization techniques in a waveform design framework 
in regard of two processing solutions, specific to OFDM radar. 

When dealing with the frequency domain processing solu¬ 
tion, we pointed out that the focus shall be the minimization 
of the PMEPR. This flgure-of-merit is a salient aspect of 
the OFDM pulse whose time domain signal is subject to 
strong variations. In this context, we suggested using a single 
objective GA. We showed that this evolutionary technique 
can produce solutions that outperform robust state-of-the-art 
methods. In particular, we challenged several limitations of 
those methods and demonstrated the merit of our approach for 
the case 1) when the spectrum is composed of non contiguous 
subcarriers and 2) when the phase codes are constrained to 
discrete values, as it is the case with specific constellations. 

Next, we explained that the time domain processing solution 
requires optimal autocorrelation sidelobe and PMEPR proper¬ 
ties. We showed that it is possible to apply another evolution¬ 
ary algorithm, namely the NSGA-II to search for optimal sets 
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of phase codes in a Pareto sense. We observed that substantial problem is equivalent to finding an optimal weight vector and 

improvement is achieved in comparison to random coding. The indicated that here again, a single GA could be used. We 

Pareto optimal solutions obtained with this technique allow characterized the enhanced detection in terms of the metric, 

various options for different design preferences. We observed namely, average power, that gives 0 dB when the OFDM pulse 

that, at a fixed time bandwidth product, the case with the has a fiat spectrum of unit energy. Then, we illustrated the 

smallest number of subcarriers produces the best Pareto front, convergence of the PMEPR and found that we could obtain 

Next, we looked at the possibility to treat the PMEPR as a optimal sets of phase codes to improve the PMEPR by about 

constrain instead, while the two objectives are the PSLR and 3 dB in comparison to the case where random codes would 

the ISLR. We demonstrated that our modified NSGA-II could be used. 

produce Pareto fronts with equivalent characteristics in terms Overall, our optimization methods run from few 10s of 
of the ISLR as for the unconstrained case. The same does not seconds to few 100s of seconds, hence real time opera- 
hold true for the PSLR, since the improvement becomes less tions would require platforms with high processing power, 
in the constrained case. Alternatively, they can be run offline. The results are then 

Ultimately, we presented a case study where we incor- stored in lookup tables for future use. Improvement of the 

porated the PMEPR optimization technique into a matched execution time for our GA based techniques and optimization 

illumination strategy. We showed that both optimization prob- of the convergence are directions for future work. Alternative 

lems could be decoupled, if the weights are reserved for the optimization strategies, among other, convex methods are also 

detection enhancement part while the phase codes are saved for of the next research directions, 

the PMEPR part. We demonstrated that the enhanced detection 




